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I.  Introduction 

•hat  is  the  flow  induced  hy  an  isolated  steady  torn*  or 
sink  benaatn  a  fro#  surface?  This  sieple  quntioa  does  not  appear 
to  hm  •  aiapla  answer.  If  tki  source  is  e  line  source  of  strength 
a,  in  two-diaensioaai  trrot at  tonal  flow  of  an  incoopressihle  laviscid 
fluid  of  infinite  depth,  and  is  situated  at  lUbaargcnce  h  beneath 
the  taidisturbed  level  of  the  free  surface  under  gravity  g,  then 
there  is  only  one  diaansiontess  pa rear ter,  the  (squared)  Pvoude 


M-o'/C**)  Cl.l) 

and  we  eight  aspect  to  find  a  solution  for  every  value  of  f* » 

In  fact .  the  problea  cannot  be  solved  without  further 
specification  of  the  nature  of  the  free ‘Surface  disturbance  taeediately 
above  the  source.  Soar  previous  investigator*  |2|,  (1)  have  assuaad 
a  stagnation  point,  and  have  sought  results  for  snail  values  of  f* , 
further  studies  of  this  type  of  flow  have  been  and#  recently  by  the 
present  authors,  and  will  he  reported  elsewhere.  A  feature  of  these 
stagnation-point  solutions  is  the  presence  of  short  waves,  which 
steepen  as  f*  increases,  end  these  eolations  soea  to  be  confined 
to  f*  <  «. 
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If  ••  fort*  (?-•>  w  keH  ii  mm  litcnt*  sit  of  values  ®  r  ®k« 
milk  #i  >  0  9*4  <  i,  ud  in  addition  require 

•otfc  Cf  lf)  wd  (2- Id)  is  held,  there  results  a  set  of  M*2  non- 
lieear  algebraic  equation*  involving  the  «*t  coefficients  ^  , 

»*•.!. 2.  -.h,  various  nueerical  net hods  can  be  used  to  solve  this 

set  of  equation*.  m  first  mm  not  decide  just  how  may  of  the 
coefficients  are  to  be  considered  as  unknown. 

If  i he  Hondo  nawber  f  is  prescribed,  then  (2,6)  deteroines 
the  leading  coefficient  b*  la  principle.  it  is  then  possible  to 
»reo»  »\  **  •  *w  *f  *  unknowns.  However,  all  atteopts 

to  solve  this  stolen  with  input  be  failed.  There  was  tone  indication 
t»ot  with  I  t  *  »e  *  If.  success  was  near,  and  it  was  then  suspected 
isos  *  «el«e  ion  mgsi  esist  onl*  for  none  oporto!  Front*  nuober. 
corses  peed  ing  to  a  h«  *otue  to  this  range .  Therefore,  the  nueericel 
.» sen oiwse*  wese  nod! fled  to  allow  b»  to  be  an  unknown  rather  than 
co  leper  O»oe« it ».  and  the  result  nos  inoediete  and  cooplete  success, 
wtth  sapid  cower  grace  to  a  solwt  ton  at  b«  *  IM.'S'.  i  e.  f*  •  12.622. 

He  oc  foci  eel  hod  it  a  vmcs  iteration,  ia  which,  if  b 

I*  os  appro  i  toot  tee  to  the  desired  solution,  then  a  better  approx  inat  ion 
it  h  «  to  ,  where  #*  it  ehtained  iw  solving 
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h 


(2.1) 


If  we  choose  *kl,  ()  II  tea  he  solved  subject  to  ;J»e  linear 
eonatsaiW's  f;  Mi  bf  oof  standard  linear-equation  package. 

It  is  set  difficult  to  obtoia  on  evpMcit  fbruula  fer  »he  Mtrlx 


element  3P/3b^  by  differentiation  of  (2.8)-(2.11) .  Uniformly 
spaced  6fc  were  found  to  be  satisfactory. 

The  iteration  process  can  be  started  with  guessed  values  such 
as  b«*1.8,  b}*0.5,  bj«0.6,  b3*-0.1,  and  all  other  coefficients  zero. 
In  practice  it  was  found  convenient  to  start  with  a  low  value  (say  5) 
for  N,  and,  once  the  iteration  converged  at  that  N,  to  use  the 
resulting  coefficients  as  a  starting  guess  for  iterations  at  a  higher 
value  of  N.  Convergence  is  very  rapid  at  any  fixed  N,  no  more 
than  S  iterations  being  ever  needed  to  reduce  the  maximum  value  of 
rCOfc  )  below  10' 5 . 

Table  1  shows  values  of  bo  from  a  run  in  which  N  was 
successively  increased  in  steps  of  5.  The  final  value  bo*l. 84257 
at  N»25  is  accurate  to  at  least  5  figures.  Table  2  shows  the 
coefficients  b  ,  j»0,l,2, . . .,10.  The  free  surface  is  shown  in 

Figure  1.  All  computations  were  carried  out  on  a  TRS-80  micro¬ 
computer. 
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TABLE  1 


N 

bo 

5 

1 .86935 

10 

1.84223 

IS 

1.84260 

20 

1.84256 

25 

1.84257 

TABLE  2 

j 

b5 

0 

1.84257 

1 

0.41325 

2 

0.55982 

3 

-0.06731 

4 

0.01766 

5 

-0.00613 

6 

0.00248 

7 

-0.00111 

8 

0.00053 

9 

-0.00027 

10 

0.00014 

[ 
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4. 


Conclusion 


We  have  provided  here  both  negative  (failure  to  achieve  solutions 
at  general  input  bo  values)  and  positive  (success  to  high  accuracy 
with  bo  as  an  unknown)  nuaerical  evidence  that  a  cusp-like  flow 
exists  only  for  a  unique  Froude  number,  close  to  F2  *  12.622. 

Several  questions  are  raised  by  this  conclusion.  If  this  cusped 
solution  exists  only  at  F2  *  12.622..,  what  happens  at  F2  »  12  or 
F2  *  13?  The  present  conclusion  relates  only  to  existence  of  a  steady 
flow, and  an  obvious  but  hardly  satisfying  answer  to  the  above  question 
is  that,  if  the  source-like  flow  is  started  from  rest,  a  steady  state 
cannot  be  achieved  if  F2  12.622...  But  then,  what  happens  instead 
of  a  steady  state? 

It  is  also  notable  that  our  present  steady  results  are 
independent  of  the  sign  of  m,  i.e.  are  as  valid  for  sinks  as  for 
sources.  In  practice,  it  seems  rather  more  likely  that  a  steady 
flow  with  a  downward -cusped  free  surface  could  occur  for  a  sink, 
than  for  a  source.  Although  the  steady  equations  are  independent  of 
the  sign  of  m,  this  is  not  so  for  the  corresponding  unsteady 
equations,  and  it  is  possible  that  the  present  solutions  could  be 
stable  for  m  <  0  but  unstable  for  m  >  0. 
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